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We consider a waveguide formed in a clean graphene monolayer by a spatially inhomogeneous 
magnetic field. The single-particle dispersion relation for this waveguide exhibits a zero-energy 
Landau-like flat band, while finite-energy bands have dispersion and correspond, in particular, to 
snake orbits. For zero-mode states, all matrix elements of the current operator vanish, and a finite 
conductance can only be caused by virtual transitions to finite-energy bands. We show that Coulomb 
interactions generate such processes. In stark contrast to finite-energy bands, the conductance is not 
quantized and shows a characteristic dependence on the zero-mode filling. Transport experiments 
thereby offer a novel and highly sensitive probe of electron-electron interactions in clean graphene 
samples. We argue that this interaction-driven zero-mode conductor may also appear in other 
physical settings and is not captured by the conventional Tomonaga-Luttinger liquid description. 

PACS numbers: 72.10.-d, 72.80.Vp, 71.10.Pm 


I. INTRODUCTION 

Electronic phases exhibiting flat bands have attracted 
considerable attention over the past few decades m- 
For instance, flat bands can arise from interference ef¬ 
fects on a geometrically frustrated lattice. On the non¬ 
interacting level, due to the lack of dispersion, one ex¬ 
pects insulating behavior when the Fermi level resides 
inside the flat band, such that the conductance vanishes 
identically at zero temperature. For lattice models host¬ 
ing almost flat bands, it is well known that interactions 
can cause dramatic effects such as topologically nontriv¬ 
ial fractional Chern insulator phases gHS]. Even topo¬ 
logically trivial phases without any dispersion can show 
remarkable behavior. For instance, in the case of long- 
range unscreened interactions, the conductance of the so- 
called Ta lattice can be finite and exhibits a highly non¬ 
trivial dependence on the filling factor [7]. Somewhat 
related conclusions have been obtained for interacting 
fermions with short-range interactions on lattices with 
geometrically frustrated unit cells, where at certain fill¬ 
ing factors the noninteracting theory predicts insulating 
behavior but repulsive interactions cause the existence of 
delocalized two-particle states Such effects have 

been studied in detail for diamond chains, where flat 
band formation emerges due to Aharonov-Bohm caging. 
However, two-particle delocalization does not necessarily 
imply that the many-particle electron system will have 
finite conductance at finite density ElE]. Similar issues 
have also been discussed for interacting bosons [14] and 
for cold-atom systems Ba¬ 
in this work, we show that a finite conductance is 
generated by Coulomb interactions in another flat-band 
system, referred to as “magnetic graphene waveguide” 
(MGW) in what follows. Our analysis for the MGW 
reveals that interactions can turn a noninteracting in¬ 
sulator into a conductor, even though electron-electron 



Figure 1. Schematic sketch of the MGW setup viewed from 
above. A magnetic field +Bez is applied everywhere in the 
graphene plane except for the waveguide region |a:| < d/2, 
where the field is reversed. Near the field switching lines, 
quantum modes corresponding to classical snake orbits are 
realized, while only Landau-quantized cyclotron orbits are 
present far away from these lines. Several classical orbits 
corresponding to n = 1 are schematically illustrated. The 
system dimensions are and Ly, respectively. 

interactions usually suppress the conductance [TMTB] . 
This effectively one-dimensional (ID) zero-mode conduc¬ 
tor falls outside the conventional Tomonaga-Luttinger 
liquid (TLL) description of interacting ID conductors 
BSj- We expect that such a state also appears in other 
physical settings and provide an in-depth description of 
its properties in a MGW. 

Our MGW setup is illustrated in Fig. where a 
clean graphene monolayer is exposed to a static inho¬ 
mogeneous magnetic field. Very long mean free paths 
have already been realized in graphene, e.g., by using 
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boron nitride as substrate |20] . Low-energy quasipar¬ 
ticles close to the neutrality point then correspond to 
massless Dirac fermions in two spatial dimensions (2D) 
[mug. The magnetic field is taken spatially inhomoge¬ 
neous along, say, the x-direction, B — B{x)ez, where we 
focus on orbital fields such that only the perpendicular 
( 2 ;-)component matters. To be specific, we study the field 
profile, see Fig. 

B{x) = B sgn(|x| — d/2) (1) 


describing a MGW of transverse width d where the mag¬ 
netic field is reversed {B —>■ —B) in the waveguide region 
dehned by |a;| < d/2. Such inhomogeneous magnetic field 
profiles allow one to guide Dirac fermions [231 - 128] . and 
single-particle spectra of the resulting MGW have been 
analyzed in detail |2SH33]- For d —>■ 0, Eq. Q reduces 
to the familiar homogeneous field case and one recovers 
the well-known dispersionless relativistic Landau levels, 
including a zero mode |34H37j . Importantly, such a zero- 
energy band (with band index n = 0) is also present for 
finite d, forming the flat band of interest below, while 
all other (n ^ 0) bands acquire dispersion. In particu¬ 
lar, near x = idj2, pairs of counterpropagating “snake 
states” develop [211 1551110] , which are either of elec¬ 

tron (n > 0) or hole (n < 0) type. Glassical snake orbits 
forming near zero-field lines are illustrated schematically 
in Fig. We note that snake states have already been 
observed experimentally in graphene |41II43] . including 
studies of the ballistic (disorder-free) limit [43] . 

When the n = 0 level is partially filled, with all 
negative-energy states occupied, the noninteracting con¬ 
ductance of the MGW vanishes identically since all n = 0 
current matrix elements are zero. We here demonstrate 
that, nonetheless, a finite conductance emerges due to 
inter-band Goulomb interactions (corresponding to Lan¬ 
dau level mixing for d = 0 [44]). The transport features 
related to the zero-energy band differ qualitatively from 
those of finite-energy bands, and thus should be easy to 
identify experimentally. 

For the guidance of the focused reader, we now sum¬ 
marize the content of the following sections before start¬ 
ing with the detailed description of our results. We 
describe the solution of the single-particle problem in 
Sec. 113 where in Landau gauge the problem is homo¬ 
geneous along the y-direction and eigenstates for given 
band index n are classified by the momentum ky = k. 
Goulomb matrix elements are discussed in Sec. TO fol¬ 
lowed by a st udy o f the theory projected to the zero-mode 
sector in Sec. Ill Employing the Hartree-Fock (HF) ap¬ 
proximation, we find that the spatial inhomogeneity of 
the magnetic field generically leads to dispersion of the 
HF single-particle energies £k ■ The HF ground state rep¬ 
resents a filled Fermi sea, where all n = 0 states with 
|fc| < kp(v) are occupied, with the zero-mode Filing 
factor V. We determine the interaction-induced Fermi 
momentum kp^v), and show that the conductance still 
vanishes within the zero-mode description. In Sec. m 
the zero-temperature linear conductance, G, will be com¬ 


puted from the Kubo formula, taking into account inter¬ 
band interactions through a systematic perturbative ex¬ 
pansion up to second order. For the zero-energy bands, 
we find completely different transport properties when 
compared to finite-energy bands. In the latter case, 
which has been studied in Ref. [3T], G is independent 
of the band filling and quantized in units of the conduc¬ 
tance quantum Gq = e^/h. Moreover, finite-temperature 
corrections to this quantized value follow the predictions 
of TLL theory. In marked contrast to these finite-energy 
bands, the conductance found for the zero-energy case 
strongly depends on the filling and is not described by 
TLL theory. Finally, in Sec. [^ we shall put our results 
into a general context. Some calculational details can be 
found in the Appendix. Throughout this paper, we focus 
on the most interesting zero-temperature limit and often 
employ units with h = 1. 

II. MAGNETIC GRAPHENE WAVEGUIDE 
MODEL 

A. Single-particle description 

We begin by discussing the single-particle description 
of the MGW, see also Refs. [351150]. The electronic low- 
energy physics in a weakly doped graphene monolayer is 
well described by massless 2D Dirac fermions [3T]. In¬ 
cluding the vector potential A encoding a static inhomo¬ 
geneous orbital magnetic held, where only the held com¬ 
ponent perpendicular to the layer matters, B = B{x)ez, 
the single-particle Hamiltonian is given by 

Hq = vpcr ■ -I- , (2) 

where vp « 10® m/s is the Fermi velocity, cr = {ax, cry) 
contains 2x2 Pauli matrices, and V = {dx,dy). The 
Pauli matrices act in the sublattice space corresponding 
to the two carbon atoms in the basis of graphene’s honey¬ 
comb lattice. The lengthscale on which B{x) changes is 
assumed large against the lattice spacing (2.46 A), such 
that valley {K point) mixing is irrelevant. The inclusion 
of spin and/or valley degrees of freedom is left to future 
work, although this step is not expected to signihcantly 
affect any of our conclusions. As a consequence, we ne¬ 
glect the spin and valley degrees of freedom in Eq. ([^ 
and focus on a single Dirac cone. 

We consider the magnetic held prohle in Eq. Q, de¬ 
scribing a MGW of width d, with constant held B(^ = B 
everywhere except in the strip |a;| < d/2 where B ^ —B, 
see Fig. [T] Using the Landau gauge, the vector potential 
is given by 

{ X + d, X < —d/2, 

—X, |a:| < d/2, (3) 

X — d, X > d/2. 

The specihc form in Eq. 0 has been adopted in order 
to keep the discussion focused and to allow for analyti- 
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cal progress. However, other field profiles creating such 
a MGW are expected to yield similar results [30], and 
indeed it is not necessary (nor even desirable) to have 
atomically sharp changes in the magnetic field profile. 

Below we often measure lengths (energies) in units of 
the magnetic length Ib (magnetic energy Eb), 

Ib = ^Jhc/eB, Eb = Hvf/Ib- (4) 

We mention in passing that instead of a true mag¬ 
netic field, one could also employ strain-induced pseudo- 
magnetic fields [451146| . and with minor modifications, 
similar physics can also be realized by employing the sur¬ 
face states of 3D topological insulators gZHlS]. 

Periodic boundary conditions along the j/-direction 
quantize the conserved momentum k = ky along the 
MGW. Spinor eigenstates, 'in,k{x,y), solving the Dirac 
equation. 


— En^k^n,k: (5) 

with eigenenergy E^^k can thus be classified by k and 
the integer band index n. For a homogeneous field 
{d = 0), Eqs. Q and ([^ lead to the well-known 
fc-independent relativistic Landau levels [2T| I34H37] . 

^n.k°'’ = sgn(n)v^^, and n can be identified with 
the Landau level index. When d is finite, single-particle 
states in general exhibit dispersion, but bands with dif¬ 
ferent band index n remain always separated by a fi¬ 
nite gap jini 1301 ■ Moreover, even in this inhomogeneous 
case, a flat zero-energy band is present, En=o,k = 0. In 
this work, we study whether this zero-mode band can 
carry electric current when Coulomb interaction effects 
are taken into account. 

Owing to momentum conservation along the y- 
direction, we have 

4'«.fe(a:,y) = ^J=ll)n,k{x), IpnA^) = txnk{x) ) ’ 

( 6 ) 

with the normalization condition f dx{(A fc + Xn fc) ~ 1- 
As shown in App. ^ both cj)n,k and Xn,k can be cho¬ 
sen real-valued, ancime single-particle problem is solved 
by matching the respective spinor solutions in the three 
regions in Fig. at a: = ±dj2. The solutions have the 
symmetry properties 

i>n,k{A = <Xzi’-n,kiA A ^ 0), (7) 

AAA = i-^r^^A,-k{-x), 



Figure 2. Dispersion relation of the MGW for d = 21 b- Only 
the marked three bands (n = —1,0,1) will be kept later on 
in Sec. |IV C| assuming that the zero mode is partially filled. 
We note that the n = 1 and n — 2 bands are separated by 
an avoided crossing not apparent on the shown scale. Energy 
(momentum) is given in units of Eb (A)’ 0. 


are given by 


In,k-n'k'=VF J dx AnAA’^V^n',k'{x), 


( 8 ) 


where the dagger denotes transposition and complex con¬ 
jugation of spinor wave functions [49]. As a consequence 
of Eq. 0. they satisfy the symmetry relations 

d^n,k;—n,k' — In,k'\ — n^k-: (9) 

A,k;n^k' — ( 1) ^Q, — k;n, — k': 


with arbitrary integer n. 

We next describe the spectrum and the eigenstates 
for the MGW, where the dispersion relation is shown in 
Fig-i for d = 21 b- As described in App. [^ the spec¬ 
trum and the eigenstates follow by numerical solution of 
an eigenvalue problem. However, analytical progress is 
possible for the zero modes. Indeed, as dictated by in¬ 
dex theorems [36l|37], there must be zero-energy states, 
En=o.k = 0, also for finite d, cf. Fig. Using units with 
Zb = 1, their analytical form is given oy 


AAA 


xoAA 


0 

Ao.kix) 

gd|x|-(a;-|-fc)^/2 


N, 


0,k 


( 10 ) 

1, \x\ > d/2, 

^{\x\-d/ 2 f^ |a;|<d/2. 


where indicates complex conjugation of both spinor 
entries. The first relation says that a particle-hole trans¬ 
formation connects states with n and —n (for n 0), and 
thus the single-particle spectrum is mirror-symmetric 
around zero energy, E-n,k = —En,k- The second relation 
follows from inversion symmetry together with the node 
rule in ID. Moreover, the matrix elements of the current 
/ (evaluated at y = 0) along the waveguide (?/-)direction 


with normalization constant A^o.fc- Due to the absence of 
an upper spinor component, all zero-mode current matrix 
elements vanish identically, 

Io,k;0,k' = 0. (11) 

For d = 0, the Xo,kix) functions reduce to shifted har¬ 
monic oscillator ground-state wavefunctions describing 
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Figure 3. Main panel: Single-particle gap Eg vs waveguide 
width d in units as in Fig. Inset: Long-wavelength part of 
the dispersion relation for n = 1 and several values of d. 


the n = 0 Landau level. For d ^ 0, the probability den¬ 
sity distribution, , has two local maxima near 

(for |A:| ^ d/l'^) but never inside the waveguide region. 
For later use, we also mention that a local minimum of 
l' 0 o,fc(^)P exists at X = kl^ when \k\ < d/2l‘^. 

Since the n ^ 0 (and particularly the n = ±1) snake 
states exhibit their probability density maximum near 
the null lines of the magnetic field [SD] at x = ±c?/ 2 , 
cf. Fig. interaction-induced transitions between n = 0 
and n ^ 0 states are therefore only important for |fc | < 


d/lg. The zero-mode conductance G discussed in Sec.lIV 
is caused by precisely such transitions. In fact, as long 
as virtual band transitions to n 7 ^ 0 states are excluded, 
G = 0 holds on general grounds since interactions cannot 
generate an upper spinor component in Eq. (10) from the 
zero-mode sector only, cf. Eq. 

Our perturbative approach in Sec. m holds as long as 
the single-particle gap, Eg = En=i,k= 0 i which separates 
the zero mode from the n = 1 band, is large against 
the typical Coulomb energy scale (see Sec. IIB). Figure 
|3 shows Eg as a function of the MGW width d. The 
^served decrease suggests to subsequently consider only 
d<2lB- Nonetheless, as shown in App.[^ it is also in¬ 
structive to discuss the limit d —?► 00 , where a rapid (ap¬ 
proximately exponential) decrease of Eg is seen in Fig. 


B. Coulomb interactions 

Next we turn to a second-quantized description and 
include Coulomb interaction effects. The fermion anni¬ 
hilation field operator at r = (x, y) is written as 

=^^ri,k(r) Cri,k, ( 12 ) 

n,k 


with fermion operators subject to the standard anti¬ 
commutator algebra, {c„ = Snn'dkk' and so on. 


and 'i’n,k in Eqs. ([^ and ([§. Using the units in Eq. 
the second-quantized interaction Hamiltonian is given by 

7L/ = - J drdr'U(r — r')^^(r)'F^(r')4'(r')'F(r), (13) 
where we consider a gate-screened Coulomb potential. 


V{r) = a I ^ — 




4i?2 


(14) 


Here a denotes graphene’s effective fine structure con¬ 
stant, with typical values in the range a « 0.1 to 2 , de¬ 
pending on the dielectric properties of the surroundings 
[m 122 ]. The second (image charge) term in Eq. ( [l4| 
comes from a parallel metallic plate, e.g., due to a gate 
electrode, at distance R from the graphene layer. The 
Fourier transform of Eq. (14) is given by (we remind the 
reader that k = ky) 


V{kx, k) = 27rQ; 


1 _ g-2flv^fcS-7fc2 


(15) 


Inserting Eq. (12) into Eq. (13) and exploiting momen¬ 


tum conservation along the y-direction, we obtain 

TT ^ •,rini,n2,n3,ni) 

ni.ra2.n3,"4 k,k' ,q 
^ ^ni,k^n2,k'^n3,k' — q'^n4,k+qi 

where the Coulomb matrix elements 


(16) 


U, 


({".}) 


1 


dk 




(ni,n4) 


ikx)4':-^;\-k.) 

(17) 


are expressed in terms of form factors. 


T' 

-J I 


{n,n') 


k^q 


(kx) = j 




(X). (18) 


We show in App. that the Coulomb matrix elements 
in Eq. 0 are real-valued and subject to the symmetry 
relations 


Tr(ni,n2,n3,"4) _ x 4(«2 ,»! ,"4 ," 3 ) 
*k,k'\q — ^k',k\ — q 

(n4,n3,n2,ni) 


(19) 


^k-\-q,k^ — q\ — q 

_ /_ ^\ni-\-n2+n^-\-n4Y-ini,n2,n3,n4) 

X / — k ^ — A) ^ j — q 

In order to obtain numerical values for the Coulomb 
matrix elements, we first compute the form factors in 
Eq. (18 1 by numerical integration over x, taking into ac¬ 


count their symmetry properties. Given the form factors, 
the remaining fca;-integration in Eq. ( |17[ ) can then be eval¬ 
uated numerically in an efficient manner. 

Notably, the zero-mode {n = n' = 0) form factors can 
be evaluated analytically, 




Y 


3/Y{k)Y {k + q) 


( 20 ) 
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Using again units with Ib = 1, we here use the complex- 
valued auxiliary function 

Y{z) = Y, (e^^^'^^'erfc (±[z T d/2]) (21) 


± ie ^^+‘'^/^erfc {-i[z =F 


with the complementary error function erfc(z) [501 151j. 
For real-valued argument z, the function Y{z) is real 


and positive. In Sec. Ill we shall also refer to the ho¬ 
mogeneous case d —)■ 0, where Eq. simplifies to 
Y(z) —>■ 2e^ . The form factors in Eq. ( [20| ) then become 

moSe Coulomb matrix elements 


', resulting in the zero- 




( 0 . 0 . 0 , 0 ) 


k^k'\q 


2a r 


d —^0 


Ly - 


Kq {\q{k -k' + g)|) 

- Ko (|<?|V4i?2 + (fc-fc' + g)2 


( 22 ) 


with the modified Bessel function Kq [50] . 

Finally, it simplifies our subsequent analysis to use an¬ 
tisymmetrized Coulomb matrix elements throughout the 
remainder of this paper, 


W, 


(ni,n2,n3."4) 


/c,/c' ;q 




■r r(ni,n 2 ,n 3 ,n 4 ) _ -r r(n 2 ,ni ,113,114 

^k,k';q ^k',k;q+k-k' 


')■ 

(23) 

which follow from Eq. ( |17| ) by antisymmetrization under 
the exchange {ni,k) o {n 2 ,k'). This antisymmetriza¬ 
tion simply reflects the fermionic anticommutator alge¬ 
bra. The matrix elements in Eq. ( p^ are also real-valued 
and enjoy the same symmetry relations, see Eq. (19), as 
the 


In the zero-mode theory, there is no kinetic energy 
term and the Hamiltonian equals Hj in Eq. (16), with all 


Uj = 0 and the form factors in Eq. (20). Unfortunately, 


numerically exact solutions for this interacting problem 
are already out of reach except for very small system 
size. Here we instead proceed by employing the textbook 
Hartree-Fock (HF) approximation [52]. However, going 
beyond this approximation is expected to cause at most 
quantitative - but not qualitative - modifications of the 
interaction-induced conductance discussed later on. Note 
also that for the corresponding homogeneous (d —>■ 0) 
problem, HF calculations give a good understanding of 
the physics away from rational filling factors related to 
the fractional quantum Hall effect [211 ESI [53li58] . As HF 
parameters we choose the occupation numbers 


nk — { Cq feCQ 


(25) 


where the expectation value is self-consistently taken 
with respect to the HF approximation of the zero-mode 
Hamiltonian. For given filling factor z/, the HF pa¬ 
rameters have to be determined under the condition 
'Thk'^k = N = vNs. By choosing the as HF parame¬ 
ters, we disregard the possibility of charge density wave 
or Wigner crystal formation [52]. Note that Wigner crys¬ 
tallization was reported in the homogeneous case (d = 0) 
with unscreened (i? —>■ oo) Coulomb interactions for cer¬ 
tain filling factors v [SSHEZl- However, for our MGW with 
externally screened interactions, see Eq. (14), we do not 
expect such phases. 

Defining single-particle energies as 


£k = > Wk^k'Uk', 


Wk,k' = 2iui™), 


(26) 


III. ZERO MODE SECTOR: HARTREE-FOCK 
THEORY 


We now consider the case of a partially filled zero 
mode, where all negative-energy bands (n < 0) are occu¬ 
pied while all positive-energy states (n > 0) are unoccu¬ 
pied. The n = 0 level has the filling factor i/ = N/Ns, 
where N particles occupy the n = 0 band and the degen¬ 
eracy degree, Ng, is given by the total magnetic flux in 
units of the flux quantum. 


Ar _ (-^^ 2d)Ly 

2nll 


(24) 


assuming a rectangular sample, see Fig. [^ The mo¬ 
mentum k = ky takes the values k = 2'Kny/Ly with 
—Ns/2 < Uy < Ns/2. This assumption of fully occupied 
(empty) bands with negative (positive) energy also holds 
for the interacting ground state as long as the typical 
Coulomb energy scale is small compared to the single¬ 
particle gap Eg, see Fig. In this section, we consider 
the zero-mode sector only and thus neglect all Coulomb 
interaction processes involving n 0 states. 


with the Coulomb matrix elements in Eq. (23), the HF 
estimate for the ground-state energy reads 




(27) 


The HF iteration starts with a normalized random initial 
distribution for ni~, where we assume Next the 


HF energies Sk in Eq. (26) are computed, where e_fc = Sk 


by virtue of the symmetry relations (19). The updated 


distribution n^, which is obtained by occupying the N 
energetically lowest states, therefore always remains even 
in k. The scheme is then iterated until convergence has 
been reached. 

We find that the HF ground-state energy converges 
quickly from above. However, there are many local en¬ 
ergy minima in occupation number space, and depend¬ 
ing on the initial configuration one may converge to 
states of widely different energy. We obtain the global 
minimum by comparing converged results for sufficiently 
many (typically a few hundred) randomly chosen initial 
states. The smooth behavior of all calculated quantities, 
such as the ground-state energy or the effective Fermi 
momentum, on the system parameters also confirms that 
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this procedure reliably finds the HF ground state. For 
the results shown here, graphene’s fine structure constant 
[see Eq. ([H] was taken as a = 0.5, with the MGW width 
d — 21 b as in Fig. To check our conclusions, we have 
performed additional calculations for other parameters, 
where the results (not shown) confirm the physical pic¬ 
ture presented in what follows. 

The zero-mode HF dispersion relation is shown in 
Fig.i For all studied filling factors, we found that the 
HF energies Sk develop a pronounced dip for small mo¬ 
menta, corresponding to states located mainly within the 
waveguide region, |A:| < d/l%. For larger jfcj, we instead 
expect an almost fiat dispersion (see below). However, 
Fig. a reveals a preferential population of states at large 
momenta, which are spatially localized near the bound¬ 
aries at x = ±Lx/2. To understand this feature, it is 
instructive to briefly study the homogeneous field case, 
d —>■ 0. Using Eqs. (p^ and (26), we observe that for 


a fc- independent distribution, Uk = v, the HE energies 
(26) are given by Sk ^ {vLyl2TT) J^l^j^dk'Wk,k'■ Ig¬ 
noring boundary effects, one could then effectively shift 
the integration variable k' to absorb fc, resulting in fc- 
independent energies Sk- However, for states localized 
near the edge of the sample, dispersion is already pre¬ 
dicted by this d = 0 result. Such “Pauli holes” near the 
sample edges are clearly observed in Eig. However, 
these boundary states play no role for the zero-mode 
conductance G, since they have no significant interac¬ 
tion matrix elements with snake states. We can there¬ 
fore safely ignore large-momentum states. In practice, we 
keep only single-particle states with jfcj < fc^, wh ere th e 
momentum cutoff is chosen as kr 


d/ll, cf. Sec. live 


Eor very small filling factor, v < 0.08, the fc = 0 mini¬ 
mum in the dispersion is above the Fermi level such that 
no small-fc modes are occupied. At larger fillings, how¬ 
ever, Fig. 1^ shows that this minimum in Sk drops below 
the Fermi level, and then evolves into a double minimum 
with increasing v. The latter feature can be understood 
from the reduced probability density \ipo^k{x)\'^ inside the 
waveguide region. This density has a local minimum at 
X = kll (for jfcj < d/2ll) and thus comes with a re¬ 
duced Coulomb repulsion cost. Finally, for v > 0.5, we 
find that the HF energies may exceed the single-particle 
gap, Sk > Eg- Since in that case our assumption of well- 
separated bands may break down, we shall focus on the 
window O.I iz 0.4 in what follows. Within this win¬ 
dow, the perturbative approach in Sec. |IV| is justified. 

The respective Fermi level then intersects £k for fc = 
i/cF with kp < d. The effective Fermi momentum, 
kpiv), and effective Fermi velocity, v(v), are depicted 
as a function of the filling factor in Fig. where v 
was obtained by linearizing £k around k = kp. It is 
worth mentioning that while the results in Figs. and 

were obtained for a = ao = 0.5, the corresponding 
results for a follow from a simple scaling argu¬ 

ment. In particular, while kp is independent of a, we 
find e[.“^ = (Qf/ao)efc°'°^ and thus = {a/. 

The converged HF results for £k can tentatively be in- 





k 



Figure 4. Self-consistent zero-mode HF energies Sk (in units 
of Eb) for various filling factors p. The red dashed horizontal 
lines denote the chemical potential at the respective filling, 
with all lower states occupied. We consider a = 0.5, d = 21b, 
R — d, Lx = 161s, and W = 400 basis states, see Eq. ( |24| . 
This parameter set corresponds to Ly = 209.4 Is, where |fe| = 
pNs/Ly = 61 represents the zone boundary. 
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V 


Figure 5. Interaction-induced Fermi momentum kp vs filling 
factor u for the parameters in Fig. The dotted curve is 
a guide to the eye only. Note the approximately linear in¬ 
dependence. The inset shows the effective Fermi velocity v vs 
ly, where v is given in units of graphene’s Fermi velocity vf- 


terpreted as signature for an interaction-induced single¬ 
particle dispersion, where the low-energy physics is gov¬ 
erned by a single pair of right- and left-movers with 
nearly linear dispersion relation and well-defined Fermi 
momentum. Nonetheless, the conductance remains zero 
unless we also include virtual band transitions to n ^ 0 
states. This statement holds true even under an exact 
treatment of interactions (beyond HF theory), since no 
upper spinor components in Eq. (101 and thus no finite 
current matrix elements can then be generated. 


IV. CONDUCTANCE 

In this section, we study the zero-temperature linear 
DC-conductance, G, of the MGW when the n = 0 band 
is partially filled. Within the HF approach in Sec. m 
intra-band n = 0 Coulomb interactions were shown to 
be responsible for an effective Fermi momentum kp = 
kp{i'), where ly is the filling factor, such that all single¬ 
particle states below the Fermi energy Ep = £kp are 
occupied. We here address the question: Are Coulomb 
interactions able to induce a finite conductance in the 
clean system? Anticipating the affirmative answer to this 
question, this feature offers a powerful novel way to di¬ 
rectly probe electron-electron interaction effects in clean 
graphene samples through transport measurements. 


A. Kubo formula 

We follow the standard Kubo linear-response formal¬ 
ism [^, expressing the linear conductance as 


G = -e^ lim 

CJ —>-0 


Imn«(a;) 


(28) 


where n^(a;) is the Fourier transform of the retarded 
current-current correlation function, 

n«(t) = -*0(t)G(t), c{t) = inmJm-m, ( 29 ) 

with i{t) = the Heaviside step function 0(t), 

and the normalized ground state | $) of the full Hamilto¬ 
nian H. In second-quantized notation, the particle cur¬ 
rent along the y-direction is described by the operator 

^ ~ ^n,k-,n',k' ^ (30) 

n,n' k,k' 

with the matrix elements in Eq. ([^. Following a se¬ 
quence of standard steps, the Fourier transform, C{uj), 
of the current-current correlator, C{t), is related to the 
imaginary part of n^(a;), and thus represents a spec¬ 
tral function for current fluctuations. Indeed, noting 
that C*{—t) = C{t) implies real-valuedness of G(w), we 
find Imn'^(a;) = —C{uj)/2. Furthermore, we note that 
G(0) = 0 because of C{—t) = —C{t). The conductance 
thus follows as 

G = 7rGo^(c. = 0), (31) 

where Go = e^//i is the conductance quantum. Writ¬ 
ing C{t) = X{t) — X{—t), we then need to evaluate the 
correlation function 


X(t) = ($|/(t)i(0)|d>). (32) 


At this point, it is useful to write the full Hamiltonian 
a.s H — Hq + W, where Hq captures not only the nonin¬ 
teracting part, cf. Eq. ([^, but also includes the n = 0 HE 
Coulomb interaction terms discussed in Sec. Ell Writing 
^0 = J2n,k ^ri,kc\ k^n,k^ the effective single-particle 
energies En^o,k = En,k and Eo.fc = Efc, all remaining in¬ 
teractions processes are then encoded by W, which in 
particular describes inter-band transitions. Eor IF = 0, 
the ground state |$) = |$o) corresponds to a Eermi sea 
with the occupation numbers 

fn,k = f{Kk), f{E) = Q{Ep - E), (33) 


where f{E) is the Fermi function taken, for simplicity, at 
zero temperature. 

In order to include IF, it is convenient to evaluate X{t) 
in Eq. (32) by using the Keldysh Creen’s function tech¬ 
nique, where the time evolution proceeds from t = —oo 
to t = -boo (forward branch, s = -b) and back from 
t = -boo to t = —oo (backward branch, s = —) |59| . 
Erom now on, we use the interaction picture, where time- 
dependent operators are denoted by I{t) = 

We then have to double all dynamical fields according to 
the branch s = ± of the Keldysh contour, I{t) —>■ Is{t) 
and so on. As a result, X{t) takes the form 


A(t) = ($ 


Tc 


^(oo)/_(t)/+(0) 




(34) 














where Tc is the time-ordering operator along the Keldysh 
contour, and the time-evolution operator reads 


S{oo) = 7c exp {—i dr sWs{^ 


s=± 


B. Diagrammatic expansion 


(35) 


(a) 


(c) 


(b) 



Our strategy will be to compute the conductance as 
perturbation series in the interaction term W, which cap¬ 
tures the effects of virtual band transitions. Going up to 
second order in W yields 


G = -h Gil) q(2) 0(^3)^ 


(36) 


Figure 6. Diagrammatic expansion of the conductance (see 
main text). Panel (a) shows the zeroth-order (m = 0) dia¬ 
gram, while (b) and (c) refer to the two first-order (m = 1) 
diagrams. 


Expanding Eq. (351 in powers of W, we obtain a corre¬ 
sponding series for X(t) = + • • •, with 

the m-th order term given by 


X(™)(t) = ^-^/ dn---dT„ 


ml 


Si • • • Sr 


X ($, 


^0 


Tc 


WsAn) ■ ■■Ws^(T^)I-(t)I+(0) $o) .(37) 


Application of Wick’s theorem to the time-ordered prod¬ 


ucts of noninteracting fermion operators in Eq. (37) al¬ 


lows one to interpret such expressions in a diagrammatic 
language. The propagator (“line”) in a given diagram 
then corresponds to the Keldysh Green’s function (with 
Keldysh indices s, s' = ±) [59] . 

- t') = -*(‘i>o|rc[c„,fc,,(f)cJ;^,,,,(t')]|<i>o), (38) 
with the Fourier-transformed components 


E - E„,fc -b is sgn {En,k - Ep^ 0+ 


rfs.-s) 


(E) = 2TTisfisE) S (a - , 


(39) 


where f{E) is the Fermi function, see Eq. (33). All con¬ 
ductance diagrams up to order m = 2 are shown in Figs.[^ 
and|^ They are constructed from the following rules: 

• Each diagram must contain two external (two- 
point) vertices representing the current operators 
I-{t) and 1+ (0). These vertices carry the respec¬ 
tive Keldysh index s = ±, and are denoted by filled 
circles in Figs. and 

• For an m-th order diagram, there are m internal 
(four-point) vertices representing the interaction 
Hamiltonian, i.e., factors sWs{t). These vertices 
are denoted by filled squares. 

• Using the Keldysh normalization condition 
($o|5'(oo)|<i)o) = 1 as well as ($|/±(t)|$) = 0, only 
connected diagrams have to be taken into account. 



Figure 7. Same as Fig. but for the second-order (m = 2) 
diagrams. 


We mention in passing that in evaluating low-order dia¬ 
grams it is useful to exploit a “selection rule” that allows 
one to discard certain contributions without detailed cal¬ 
culation. Indeed, owing to the vanishing current matrix 
elements in Eq. (11), a diagrammatic contribution is zero 


whenever a current vertex with Keldysh index s = ± is 
connected to two other vertices that both have the oppo¬ 
site Keldysh index —s jSD]. 

At this stage, it is convenient to introduce several aux¬ 
iliary matrices. First, we define the current matrix I 
through its matrix elements In,k-,n',k' in Eq- Second, 
the hybridization matrix A with 


A. 


n,fc;n' 


En.k 


'n' ,k 


n \ ,/li 


Jni,k 


-rxr(n,ni,ni,n') 

^^k,ki;q^0 


(40) 

describes interaction-induced virtual transitions between 
different bands. Here the intermediate summation comes 
from a fermion loop. We notice that A is diagonal in k- 
space and hence can indeed be interpreted as hybridiza¬ 
tion matrix. The commutator of the above matrices is 
given by 


57= [/,A]_, 


(41) 
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which captures the current matrix renormalization by 
virtual interband transitions to leading order. Finally, 
we introduce a fluctuation matrix K, with matrix ele¬ 
ments 


Kr, 


= E 


^k2,ki-\-k' —k 


fni,ki _ fn2,k2 


ni,ki;n2,k2 




'Ei2-,k2 


-E, 


'n\ ,fci 


V T ^j,An2,n,ni,n') 

^ J^ni,ki-,n2,k2^^ k2,k-k-ki ' 


(42) 


Below, the “o” symbol denotes Fermi level convolution, 
i.e., the matrix elements oi Ao B are given by 

{A O B'ji^^k^n',k' ~ ^ ) ^n,fc;0,ibfcF 


and “Tr^?” denotes a Fermi level trace, 

TvpA = ^ ^0,±feF;0,±fcF ■ 


(44) 


Let us now start with the zeroth-order conductance 
contribution, where inter-band transitions are absent. 
From the discussion in Sec. m this term is expected 
to vanish identically, = 0. There is only a single 
TO = 0 diagram, represented by the polarization bubble 
inFig.|^a). By virtue of Eqs. (37) and (38), this diagram 
leads to the correlation function 

X(°)(t)= ^ (45) 


Performing a Fourier transformation, we find 
l(°)(cc) = 27r ^ 

n,k;n' ,k' 

^ (1 fn,k^fn',k'^i^n,k ^n',k' 


(46) 


and with Eq. (31), the respective zero-temperature con- 


worth mentioning that for a conventional system with 
non-zero current matrix elements, e.g., the finite-energy 
bands for our MGW, Eqs. (|47|) and (48| yield finite re¬ 


sults, corresponding to the “noninteracting” conductance 
and the ballistic version of the “interaction correction” 
respectively. 

In order to encounter a finite conductance in our zero¬ 
mode system, we have to go up to second order (to = 2). 
All topologically distinct to = 2 diagrams are shown in 
Fig. [3 Using the above selection rule, we find that dia¬ 
gram (c) also gives no conductance contribution. More¬ 
over, diagrams (f) and (g) vanish as well since they in¬ 
volve products of more than two Fermi factors which 
never satisfy the resulting energy constraints. Within 
second-order perturbation theory, the conductance is 
thus obtained from the remaining diagrams in Fig. [3 
G = G^^) = Gq -I- Gb -I- Gd -I- Ge, which yield a finite 
result. Evaluating diagrams (a) and (b) together gives 


(G, + Gb)/Go = 27t^p^o TvFiSIoSI). 

Similarly, diagram (d) yields 

Gd/Go = -47rVo TiFiSIoK), 
while diagram (e) produces the contribution 
GelGo = 2TT^pl TrFiKoK). 


(49) 


(50) 


(51) 


Collecting all diagrams, we obtain a manifestly positive 
and general result for the zero-mode conductance, 

G/Go = 27rVo TrF [(<5/ - K) o {SI - K)] . (52) 


ductance contribution is 


We note in passing that Eq. (52) does not apply for d —>■ 
0, where the zero-mode dispersion Sk becomes flat and 
hence kF is not defined anymore. 


C. Zero-mode MGW conductance 


G(o)/Go = ^ Tr;^(/o/). 


(47) 


Here po = 1 /ttv denotes the density of states, with 
V = dkSk=kF, and the symbols “o” and “Tri?” were de¬ 
fined in Eqs. ([d^) and (44), respectively. Since all current 


matrix elements appearing in Eq. (47) vanish by virtue 
of Eq. (ITll) , this calculation confirms the expected result 

G(o) = a 

Let us then turn to the first-order (to = 1) diagrams. 
After some algebra, the corresponding conductance con¬ 
tribution takes the form 

GF)/Go = -27r^ pI Tr^ [/ o {SI - K)] , (48) 

where the (5/-term comes from diagram (b) and the K- 
term from diagram (c) in Fig. respectively, see also 
Eqs. (41) and (42). Because of the vanishing zero-mode 
current matrix elements in Eq. ( |11[ ), also this conduc¬ 
tance contribution vanishes, G*-^^ = 0. However, it is 


We next discuss the zero-mode MGW conductance pre- 


dicted by Eq. (52), adopting the parameters in Sec. HI 


Our main goals are (i) to reliably demonstrate the ex¬ 
istence of a finite zero-mode conductance, and (ii) to 
clarify its filling dependence. In order to simplify the 
numerical evaluation, which is quite cumbersome due to 
the presence of interaction matrix elements connecting 
all different bands, we shall here evaluate Eq. (52) by 


taking into account only the three bands n = —1,0,1 
sketched in Fig. Indeed, the n = ±1 bands are ener¬ 
getically closest to the n = 0 modes and therefore pro¬ 
duce the main conductance contribution. Moreover, to 
avoid spurious finite-size effects, we introduce a momen¬ 
tum bandwidth kc restricting the single-particle Hilbert 
space to states with |fc| < k. 
below, where 
k^ = 
range 1.5 
results. 


see Sec. HI For the results 


d = 21b, we chose the momentum cutoff 
However, taking other values within the 
^Jb ^ l.f also gave essentially identical 
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Figure 8 . Zero-mode waveguide conductance G (in units of 
Go = e^lh) vs filling u, see Eq. (521, for two values of the 
fine structure constant a. Dotted curves are guides to the eye 
only. All other parameters are as in Fig. The inset shows 
the a-dependence at two selected fillings. 


The resulting zero-temperature conductance is illus¬ 
trated in Fig. The main panel shows that the conduc¬ 
tance strongly depends on the zero-mode filling factor v, 
with a pronounced minimum around v = I'min, where we 
find r'min ~ 0.145 at a = 0.5. Near this minimum, G be¬ 
comes very small. The existence of the minimum in G{i') 
can be rationalized by analyzing the Coulomb-assisted 
hybridization of zero modes with the n = ±1 bands. For 
small filling v, and therefore small chemical potential, see 
Sec. |IIH the band n = — 1 still remains close to the n = 0 
states. In that case, the current renormalization effects 
encoded in SI are found to dominate the conductance in 
Eq. ( [5^ . With increasing filling, this conductance contri¬ 
bution begins to weaken, while the Fermi level gets closer 
and closer to the n = 1 band. Eventually, at large filling, 
the fluctuation matrix K instead dominates in Eq. (52|. 
We then encounter a nearly perfect cancellation of the4/ 
and K terms for filling factor iz « v^nin- 


The inset of Fig. shows the a-dependence of the con¬ 
ductance for two fillings v, chosen below and above i'min, 
respectively. Two comments are in order here. First, 
the limit a —>■ 0 seems to result in a finite conductance. 
This may come as a surprise, since in the absence of 
interactions, we know that G = 0 must hold. How¬ 
ever, one should keep in mind that our HF approach in 
Sec. |Tn] applies only if the frequency scale oj at which the 
conductance is probed remains well below the effective 
bandwidth of the zero mode. Since this Sk bandwidth is 
proportional to a, it vanishes in the limit a —> 0 , and 
the w —>■ 0 limit implicit in Eq. ( |52[ ) cannot be taken 
anymore. Second, the different a-dependence at the two 
fillings observed in the inset of Fig. [^simply reflects the 
fact that r'min increases with increasing a, cf. the main 
panel. This shift of t'min in turn can be rationalized by 
noting that the chemical potential also moves up with in¬ 
creasing a, and thus the cancellation point between the 


SI and K contributions is slightly shifted towards larger 
fillings. 

It is instructive to contrast the above results for 
the zero-mode conductance with the conductance found 
when the chemical potential intersects one of the n yf 0 
bands instead. This case has been studied in detail in 
Ref. [3T] , where a completely different behavior has been 
reported. When contacted by wide electrodes, the zero- 
temperature conductance then assumes a quantized value 
(in units of Go), where transport proceeds predominantly 
through snake states. This quantization can be rational¬ 
ized by noting that right- and left-moving snake states 
are spatially separated [3T]. In addition, no pronounced 
dependence of the conductance on the respective band 
filling is expected, in marked contrast to the zero-mode 
case shown in Fig. As we elaborate further in Sec. [Vj 
the finite-energy bands correspond to a realization of the 
conventional TLL phase, which could also be detected 
through the predicted power-law corrections to the con¬ 
ductance at finite temperatures m- 

To summarize, the strong dependence of the zero-mode 
conductance G{iy) on the filling factor u, see Fig. ^ 
should allow for a clear experimental signature of the 
predicted interaction-induced conducting phase. As we 
argue in the next section, such a state is distinct from a 
TLL state, and hence also not described by Fermi liquid 
theory. 


V. DISCUSSION 


In strictly ID band metals, it is well known that low- 
energy excitations near the Fermi points at ±kp are 
severely restricted as a consequence of phase-space limi¬ 
tations |52j . and therefore any nonzero electron-electron 
interaction slightly destabilizes the Fermi liquid |19| . 
The resulting phase is commonly coined “Tomonaga- 
Luttinger liquid” (TLL), where all low-energy properties 
of the system are fully determined by just two parame¬ 
ters when disregarding the spin sector [6TH63]. If Galilei 
invariance holds in addition, which is the case for the 
continuum model considered here, a single TLL param¬ 
eter, ^tll, remains. The Fermi liquid case is recovered 
for 5 tll = 1- The Kubo conductance of an infinitely 
long and clean TLL is G = gTLhGo [SIIj and single¬ 
particle correlation functions exhibit power-law behav¬ 
ior with exponents controlled by ^tll [a Eg. To give 
just one example, the equal-time single-particle Green’s 
function, cf. Eq. ( [T^ , has the asymptotic power-law de¬ 
cay (4'(a:,y)^'l'(a;,0)) -- \y\~'^ with 7 = 1 -fi [^tll + 
1/ffTLL — 2]/4. The value of 5 tll is fixed by the inter¬ 
action strength through the ground-state compressibility 

m. 


ffTLL = 


TT d'^iSojLy) 

V dkl 


- 1/2 


(53) 


where £q is the ground-state energy and v = \dkSk=kp\ 
the single-particle velocity, see also Ref. [SS]. Precisely 
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this scenario has previously been identified in our MGW 
for all n 7 ^ 0 snake states |3T]. Without interactions, 
snake states propagate uniformly at the Fermi velocity vp 
of the graphene host. They represent spatially separated 
chiral branches located near either of the two parallel 
zero lines of the magnetic field. The value of ^tll is 
then governed by Coulomb interactions between these 
oppositely moving branches and can be tuned directly 
via the MGW width d. 


On the other hand, the partially filled n = 0 MGW 
band investigated here clearly does not fit into the above 
TLL framework. First, without interactions there is no 
Fermi surface, since the zero-energy level is strictly flat. 
Second, when accounting for intra-band Coulomb inter¬ 
actions on the level of the HF approximation, we found an 
interaction-induced dispersion, where the effective single¬ 
particle energies allow us to define a Fermi momentum, 
kp, and a Fermi velocity, v, see Fig. One may then 
naively conclude that once again a TLL emerges, where 


5 TLL follows from Eqs. (27) and (53). However, this value 


of 5tll does not describe correctly the conductance of the 
system which is zero. We stress that the vanishing con¬ 
ductance of the n = 0 band alone holds true even when 
performing an exact calculation, see Sec. Ill Therefore, 
contrary to all n 7 ^ 0 bands, we conclude that the n = 0 
band cannot be described by a bosonized Gaussian field 
theory underlying the Luttinger liquid concept [S]. Only 
when accounting for virtual interband transitions to con¬ 
ducting n 7 ^ 0 bands, the n = 0 band acquires a nonzero 
conductance for which we find a quite peculiar depen¬ 
dence on the filling p. 

We expect that apart from the MGW studied here, 
similar behavior should be observable also in other set¬ 
tings. For instance, consider metallic carbon nanotubes 
with a magnetic field applied perpendicular to the tube 
axis [55]. This field should be inhomogeneous on the 
scale of the tube radius, such that a non-zero net mag¬ 
netic flux $ penetrates the tube, since only then a fi¬ 
nite degeneracy, ^/(hc/e), of the n = 0 Landau level 
is guaranteed by index theorems [35], see also Ref. [55] . 
Accounting for intra-Landau level interactions, we expect 
that dispersion of the n = 0 level is created. Nonetheless, 
the system will still exhibit insulating behavior when the 
Fermi level is close to neutrality, and only when including 
inter-Landau level interactions, a non-zero conductance 
can emerge. Similar to the case of the MGW, the actual 
value of the conductance will then give direct informa¬ 
tion on the strength of the Coulomb interaction, and the 
conductance behavior of the n = 0 level should again 
significantly differ from all n 7 ^ 0 bands. 

Our theory assumes that one works with a ballistic 
(disorder-free) sample. We expect that very weak dis¬ 
order will not qualitatively change the scenario outlined 
above, but strong disorder will introduce localization and 
thereby destroy the physics described here. Since ballistic 
transport is nowadays reachable in high-quality graphene 
devices, our results should be testable in the near future. 

To conclude, the zero-energy levels in a clean magnetic 


graphene waveguide are predicted to display qualitatively 
different conductance features than the bands of nonzero 
energy. For the latter, the zero-temperature conductance 
is quantized in multiples of the conductance quantum. 
In contrast, the zero-mode conductance is non-universal 
with a strong dependence on the filling factor. Since 
without Goulomb interactions this conductance vanishes 
identically, transport experiments offer a direct interac¬ 
tion probe. We hope that our predictions can soon be put 
to an experimental test and will inspire further studies 
of this novel transport regime. 
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Appendix A: Spectrum of MGW 


In this Appendix, we provide some details concerning 
Sec. II Aj For given fc, using the natural units in Eq. Q, 
Hq in Eq. ([^ reduces to the ID Hamiltonian 


= —iaxdx -f [fc -I- A{x)\a. 


V’ 


(Al) 


where A{x) in Eq. ® is antisymmetric under x-inversion, 
. (k) 

TZx ■ X —> —X. Due to this property, Hq ' exhibits inver¬ 
sion symmetry. 


H. 


(k) ^ 


= 0, S = TZk.TZxC, 


(A2) 


where TZ^ : k ^ —k inverts k and C denotes complex 
conjugation. The operator S has eigenvalues ^ = ±, and 
Eq. (A2) implies that the eigenstates in Eq. ^ obey the 
relation 


= -ipn-ki-x) = ±V'n,fc(a;)- 


(A3) 


Next, we note that for x < —dj2, the general solution 
at energy E = En^k can be written in the form |25j 

f D_i+e2/2 + k + d)) 

I -^^Dp 2^2 {-V2ix + k + d)) 

(A4) 

with a coefficient and the parabolic cylinder func¬ 

tion Dp{x) [SO] [51]. Similarly, with coefficients 
the solution in the waveguide region |x| < d/2 reads 

/ r \ ( Dp2f2{^V^{x—k)) \ 

/^nMIlix) - \^^l^D_,+E2/2{±V2{x-k)) J’ 

(A5) 
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while for x > d/2, one finds 


; l \ I 1 +-E 2/2 fV ^(,3^ + K ci)j 

\^-^DE2/2{V2{x + k-d)) 

(A6) 

The symmetry relations in Eq. ( |A3| n ow connect the 
coefficients in Eqs. (A4), (A5) and ( |A6[ ). With ^ = ±, 
we find the relations 


{V2{x + k — d)) 


CinM ^ ^ = e ( 


cth')'- 


(A7) 


Let us then choose the overall phase of each eigen¬ 
state such that all are real-valued. Hence 


Eq. (A71 implies that all coefficients in the vector 

c„., = are also real¬ 

valued, which in turn confirms that (/„ ^ and Xn k in 
Eq. § can indeed be chosen real-valued. 

Using the real-valued coefficient vector Cn,k, the 
matching conditions can be written in compact form as 
Mn^k{E) ■ Cn,k = 0, where the 4x4 matrix Mn^k{E) 
is easily read off from Eqs. (A4), (A5| and (|A6|). The 


eigenenergies E = En,k then follow from the condition 


det Mn^k{E) = 0, 


(AS) 


and the normalized eigenstates ipn,kix) are determined 
from the correspondi ng e igenvectors Cn,k- In general, 
the solutions to Eq. (A8) have to be obtained by using 


numerical root-finding methods ISO]. 


Appendix B: On the limit d —> oo 

Here we briefly discuss the large-d behavior of the 
single-particle solutions in Sec. |H A| In effect, increasing 
the MGW width d is a way to reverse the direction of the 
magnetic field, Be^ —?> —Be^, in the bulk of the sample. 
Eventually, as d —>■ oo, all eigenstates must approach to 
Landau levels with index n S Z of the time-reversed sys¬ 
tem. 


where the (pn are normalized eigenstates of the ID har¬ 
monic oscillator and sgn(n) = (1,0,—1) for {n > 0,n = 
0,n < 0). In particular, a zeroth Landau level must 
arise, h = 0, where now only the upper spinor component 
is nonzero as compared to Eq. (10). This development 
is nicely tracked from the inset of Fig. together with 
Fig-i Indeed, as d increases, the n = ±1 snake lev¬ 
els successively flatten in order to ultimately join at zero 
energy. (Note that the avoided crossing in Fig. shifts 
towards bigger |fc| values with increasing d.) Eventually, 
the former n = ±1 snake levels merge at zero to form the 
new zeroth Landau level, see Eq. ( |B1[ ) with tt, = 0. In 
fact, using properties of the parabolic cylinder functions 
[SnilM], one finds [see Eq. (A5|] 


i’n,ki^) 


^-l/4^-(x-kf/2 


E^O. (B2) 


Appendix C: Form factor symmetries 

In this Appendix, we study general properties of the 
form factors defined in Eq. (18). Using Eq. 0 , we find 
that they obey the symmetry relations 


4Tik.) = -4°r\k.), 


k,q 


(Cl) 




where the first two relations only hold for n ^ 0 and 
n' 7 ^ 0, while the last two are valid for arbitrary n,nL 
Furthermore, as a result of 


T 

-J I 


iT'^ikx) = {E4^'\-kx)y, (C2) 


i’h,k{x) = 


1 


a/1 -bsgn(|h|) 


ip\n\{x + k) 

isgn{n)(p\fi_i\{x -b k) 


all Coulomb matrix elements in Eq. 0 are real-valued. 
Using Eq. (Cl), we hnd that they obey the symmetry 


(HI) relations quoted in Eq. (19). 
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